New Mexico Tech 
(September 2005) 



to be published in: 

Proceedings of the Workshop "Krzysztof Woj- 
ciechowski 50 years - Analysis and Geometry 
of Boundary Value Problems" , Roskilde, Den- 
mark, 20-22 May, 2005 

Non-Laplace type Operators 
on Manifolds with Boundary 

Ivan G. Avramidi 

Department of Mathematics 
New Mexico Institute of Mining and Technology 
Socorro, NM 87801, USA 
E-mail: iavramid@nmt.edu 



Non-Laplace type operators are second-order elliptic partial differential operators act- 
ing on sections of vector bundles over a manifold with boundary with a non-scalar leading 
symbol. Such operators appear, in particular, in models of non-commutative gravity theo- 
ries, when instead of a Riemannian metric there is a matrix valued self-adjoint symmetric 
two-tensor that plays the role of a "non-commutative" metric. We construct the corre- 
sponding heat kernel and compute its first two spectral invariants. 
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1 Introduction 

Elliptic differential operators on manifolds play a very important role in math- 
ematical physics, geometric analysis, differential geometry and quantum theory. 
Of special interest are the resolvents and the spectral functions of elliptic opera- 
tors; the most important spectral functions being the trace of the heat kernel and 
the zeta function, which determine, in particular, the functional determinants of 
differential operators (see [29, 19, 18, 6] and the reviews [5, 7, 43]). 

In particular, in quantum field theory and statistical physics the resolvent deter- 
mines the Green functions, the correlation functions and the propagators of quan- 
tum fields, and the functional determinant determines the effective action and the 
partition function (see, for example [6]). In spectral geometry, one is interested, 
following Kac [36], in the question: "Does the spectrum of the scalar Laplacian 
determine the geometry of a manifold" or, more generally, "To what extent does 
the spectrum of a differential operator on a manifold determine the geometry of 
the manifold?" Of course, the answer to Kac's and other questions depends on the 
differential operator. Most of the studies in spectral geometry and spectral asymp- 
totics are restricted to so-called Laplace type operators. These are second-order 
partial differential operators acting on sections of a vector bundle with a positive 
definite scalar leading symbol. 

Since, in general, it is impossible to find the spectrum of a differential operator 
exactly, one studies the asymptotic properties of the spectrum, so-called spectral 
asymptotics, which are best described by the asymptotic expansion of the trace of 
the heat kernel. If L : C°°{V) — > C°°{V) is a self-adjoint elliptic second-order par- 
tial differential operator with a positive definite leading symbol acting on smooth 
sections of a vector bundle V over a compact n-dimensional manifold M, then 
there is a small-time asymptotic expansion as ? — > [33, 29] 

CO 

Tr i2 exp(-rL) ~ (An)-"'^ J] • (1-1) 

The coefficients At are called the global heat invariants (in mathematical literature 
they are usually called the Minakshisundaram-Pleijel coefficients; in physics liter- 
ature, they are also called HMDS (Hadamard-Minakshisundaram-De Witt-Seeley) 
coefficients, or Schwinger-De Witt coefficients). 

The heat invariants are spectral invariants of the operator L that encode the 
information about the asymptotic properties of the spectrum. They are of great 
importance in spectral geometry and find extensive applications in physics, where 
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they describe renormalization and quantum corrections to the effective action in 
quantum field theory and the thermal corrections to the high-temperature expan- 
sion in statistical physics among many other things. They describe real physical 
effects and, therefore, the knowledge of these coefficients in explicit closed form 
is important in physics. One would like to have formulas for some lower-order 
coefficients to be able to study those effects. 

The proof [33, 29] of the existence of such an asymptotic expansion (1.1) has 
been a great achievement in geometric analysis. Now it is a well known fact, at 
least in the smooth category for compact manifolds. This is not the subject of 
our interest. The main objective in the study of spectral asymptotics (in spectral 
geometry and quantum field theory) is, rather, the explicit calculation of the heat 
invariants A/, in invariant geometric terms. 

The approach of Greiner and Seeley [33, 41] is a very powerful general ana- 
lytical procedure for analyzing the structure of the asymptotic expansion based on 
the theory of pseudo-differential operators and the calculus of symbols of opera- 
tors (we will call it symbolic approach for symplicity). This approach can be used 
for calculation of the heat invariants explicitly in terms of the jets of the symbol 
of the operator; it provides an iterative procedure for such a calculation. However, 
as far as we know, because of the technical complexity and, most importantly, 
lack of the manifest covariance, such analytical tools have never been used for the 
actual calculation of the explicit form of the heat invariants in an invariant geo- 
metric form. As a matter of fact, the symbolic method has only been used to prove 
the existence of the asymptotic expansion and the general structure of the heat in- 
variants (like their dependence on the jets of the symbol of the operator) (see the 
review [7] and other articles in the same volume, and [43, 38]). To the best of our 
knowledge there is no exlicit formula even for the low-order coefficients A\ and 
A2 for a general non-Laplace type operator. 

The development of the analysis needed to discuss elliptic boundary value 
problems is beyond the scope of this paper. We shall simply use the well known 
results about the existence of the heat trace asymptotics of elliptic boundary value 
problems from the classical papers of Greiner [33] and Seeley [41] (see also the 
books [35, 19]). Our approach can be best described by Greiner's own words [33], 
pp. 165-166,: "the asymptotic expansion can be obtained by more classical meth- 
ods. Namely, one constructs the Taylor expansion for the classical parametrix [of 
the heat equation] . . . and iterates it to obtain the Green's operator. This yields, at 
least formally, the asymptotic expansion for [ the trace of the heat kernel ] ". This is 
the approach exploited in [39] for a Laplace type operator and it is this approach 
that we will use in the present paper for non-Laplace type operators. However, 
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contrary to [33, 41] we do not use any Riemannian metrics but, instead, work di- 
rectly with densities, so that our final answers are automatically invariant. Greiner 
[33], pp. 166, also points out that "Of course, at the moment it is not clear which 
representation will yield more easily to geometric interpretation." 

In spectral geometry as well as in physics the motivation and the goals of 
the study of spectral asymptotics are quite different from those in analysis. The 
analytic works are primarily interested in the existense and the type of the asymp- 
totic expansion, but not necessarily in the explicit form of the coefficients of the 
expansion. In spectral geometry one is interested in the explicit form of the spec- 
tral invariants and their relation to geometry. One considers various special cases 
when some invariant topological and geometrical constraints are imposed, say, 
on the Riemannian structure (or on the connection of a vector bundle). Some of 
these conditions are: positive (negative, or zero) scalar curvature, or positive (neg- 
ative) sectional curvature, Ricci-flat metrics, Einstein spaces, symmetric spaces, 
Kaehler manifolds etc. Such conditions lead to very specific consequences for the 
heat invariants which are obvious in the geometric invariant form but which are 
hidden in a non-invariant symbolic formula obtained in local coordinates. For ex- 
ample, if the scalar curvature is zero, then for the Laplacian on a manifold without 
boundary A2 = 0. Such a conclusion cannot be reached until one realizes that the 
integrand of A2 is precisely the scalar curvature. There are, of course, many more 
examples like this. 

Another property that does not become manifest at all in the symbolic ap- 
proach is the behavior of the heat invariants under the conformal transformation 
of the Riemannian structure and the gauge transformations. This is a very impor- 
tant property that is heavily used in the functorial approach [22, 24], but which 
is not used at all in the symbolic approach. For conformally covariant operators 
the symbolic calculus is exactly the same as for non-conformally covariant ones 
with similar results because the conformal covariance only concerns the low-order 
terms of the symbol but not its leading symbol. However, the conformal invari- 
ance leads to profound consequences for the heat invariants, zeta-function and the 
functional determinant (see [20]). 

The calculation of the explicit form of the heat invariants is a separate impor- 
tant and complicated problem that requires special calculational techniques. The 
systematic explicit calculation of heat kernel coefficients was initiated by Gilkey 
[28] (see [29, 43, 38, 7] and references therein). A review of various algorithms 
for calculation of the heat kernel coefficients is presented in [4]. The two most ef- 
fective methods that have been successfully used for the actual calculation of the 
heat invariants are: 1) the functorial method of Gilkey and Branson [22, 24, 29], 
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which is based on the invariance theory, behavior of the heat trace under con- 
formal transformations and some special case calculations, and 2) the method of 
local Taylor expansion in normal coordinates (which is essentially equivalent to 
the geometric covariant Taylor expansions of [2, 1]). The results of both of these 
methods are directly obtained in an invariant geometric form. The symbolic cal- 
culus approach, despite being a powerful analytical tool, fails to provide such 
invariant results. It gives answers in local coordinates that are not invariant and 
cannot be made invariant directly. For high-order coefficients the problem of con- 
verting such results in a geometric invariant form is hopeless. One cannot even 
decide whether a particular coefficient is zero or not. 

One of the main problems in the study of spectral asymptotics is to develop 
a procedure that respects all the invariance transformations (diffeomorphisms and 
gauge transformations in the physics language) of the differential operator. Sym- 
bolic calculus gives an answer in terms of jets of the symbol of the operator in 
some local coordinates. Thus there remains a very important problem of convert- 
ing these local expressions to global geometric invariant structures, like polyno- 
mials in curvatures and their covariant derivatives. For a general coefficient this 
problem becomes unmanagable; it is simply exponentially bad in the order of the 
heat kernel coefficient. The number of the jets of the symbol is much greater than 
the number of invariant structures of given order. This problem is so bad that it 
is, in fact, much easier to compute the coefficients by some other methods that di- 
rectly give an invariant answer than to use the results of the symbolic approach. To 
our knowledge, none of the results for the explicit form of the spectral invariants 
were obtained by using the symbolic calculus. 

Every problem in geometric analysis has two aspects: an analytical aspect and 
a geometric aspect. In the study of spectral asymptotics of differential operators 
the analytic aspect has been succesfuUy solved in the classical works of Greiner 
[33] and Seeley [41] and others [35, 19]. 

The geometric aspect of the problem for Laplace type operators is now also 
well understood due to the work of Gilkey [28] and many others (see [29, 43, 38, 
18, 7]). The leading symbol of a Laplace type operator naturally defines a Rie- 
mannian metric on the manifold, which enables one to employ powerful methods 
of differential geometry. In other words, the Riemannian structure on a manifold 
is determined by a Laplace type operator. We take this fact seriously: geometry 
(Riemannian structure) is determined by analysis (differential operator). In some 
sense, analysis is primary and geometry is secondary. What kind of geometry is 
generated does, of course, depend on the differential operator. A Laplace type 
differential operator generates the Riemannian geometry. 
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As a result, much is known about the spectral asymptotics of Laplace type 
operators, both on manifolds without boundary and on manifolds with boundary, 
with various boundary conditions, such as Dirichlet, Neumann, Robin, mixed, 
oblique, Zaremba etc. On manifolds without boundary all odd coefficients vanish, 
Mk+i = 0, and all even coefficients A2k up to Ag have been computed in our PhD 
thesis [1], which was published later as a book [6] (see also [29, 2, 16, 11, 44], the 
reviews [5, 7, 43] and references therein). By using our method [2] Yajima et al. 
[44] computed the coefficient Aio recently. Of course, this remarkable progress 
can only be achieved by employing modern computer algorithms (the authors of 
[44] used a Mathematica package). The main reason for this progress is that the 
heat kernel coefficients are polynomial in the jets of the symbol of the operator 
(which can be expressed in terms of curvatures and their covariant derivatives); it 
is essentially an algebraic problem. 

On manifolds with boundary, the heat invariants depend on the boundary con- 
ditions. For the classical boundary conditions, like Dirichlet, Neumann, Robin, 
and mixed combination thereof on vector bundles, the coefficients have been 
explicitly computed up to A5 (see, for example, [37, 22, 24, 3]). 

A more general scheme, called oblique boundary value problem [34, 32, 31], 
which includes tangential derivatives along the boundary, was studied in [15, 17, 
16, 26, 27]. This problem is not automatically elliptic like the classical bound- 
ary problems; there is a certain condition on the leading symbol of the boundary 
operator that ensures the strong ellipticity of the problem. As a result, the heat 
invariants are no longer polynomial in the jets of the boundary operator, which 
makes this problem much more difficult to handle. So far, in the general case only 
the coefficient Ai is known [15, 17]. In a particular Abelian case the coefficients 
A2 and A3 have been computed in [27]. 

A discontinuous boundary value problem, the so-called Zaremba problem, 
which includes Dirichlet boundary conditions on one part of the boundary and 
Neumann boundary conditions on another part of the boundary, was studied re- 
cently in [1 1, 42, 25]. Because this problem is not smooth, the analysis becomes 
much more subtle (see [11, 42] and references therein). In particular, there is a 
singular subset of codimension 2 on which the boundary operator is discontinu- 
ous, and, one has to put an additional boundary condition that fixes the behavior at 
that set. Seeley [42] showed that there are no logarithmic terms in the asymptotic 
expansion of the trace of the heat kernel, which are possible on general grounds, 
and that the heat invariants do depend on the boundary condition at the singular 
set; the neglect of that simple fact lead to some controversy on the coefficient A2 
in the past until this question was finally settled in [42, 11]. 
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Contrary to the Laplace type operators, there are no systematic effective meth- 
ods for an explicit calculation of the heat invariants for second-order operators 
which are not of Laplace type. Such operators appear in so-called matrix geom- 
etry [8, 9, 10, 12], when instead of a single Riemannian metric there is a matrix- 
valued symmetric 2-tensor, which we call a "non-commutative metric". Matrix 
geometry is motivated by the relativistic interpretation of gauge theories and is 
intimately related to Finsler geometry (rather a collection of Finsler geometries) 
(see [8, 9, 10]). For an introduction to Finsler geometry see [40]. 

Of course, the existence and the form of the asymptotic expansion of the heat 
kernel is well established for a very large class of operators, including all self- 
adjoint elliptic partial differential operators with positive definite leading symbol; 
it is essentially the same for all second-order operators, whether of Laplace type 
or not. However, a non-Laplace type operator does not induce a unique Rieman- 
nian metric on the manifold. Of course, one can pick any Riemannian metric 
and work with it, but this is not natural; it does not reflect the properties of the 
difi'erential operator and its leading symbol. Therefore, it is useless to try to use 
a Riemannian structure to cast the heat invariants in an invariant form. Rather, 
a non-Laplace type operator defines a collection of Finsler geometries (a matrix 
geometry in the terminology of [8, 9, 10, 12]). Therefore, it is the matrix geome- 
try that should be used to study the geometric structure of the spectral invariants 
of non-Laplace type operators. This fact complicates the calculation of spectral 
asymptotics significantly. Of course, the general classical algorithms described in 
[33,41] still apply. 

Three decades ago Greiner [33], p. 164, indicated that "the problem of inter- 
preting these coefficients geometrically remains open". There has not been much 
progress in this direction. In this sense, the study of geometric aspects of spectral 
asymptotics of non-Laplace type operators is just beginning and the corresponding 
methodology is still underdeveloped in comparison with the Laplace type theory. 
The only exception to this is the case of exterior p-forms, which is pretty simple 
and, therefore, is well understood now [30, 23, 21]. Thus, the geometric aspect of 
the spectral asymptotics of non-Laplace type operators remains an open problem. 

A first step in this direction was made in our papers [13, 14]. We studied a sub- 
class of so-called natural non-Laplace type operators on Riemannian manifolds, 
which appear, for example, in the study of spin-tensor quantum gauge fields. The 
natural non-Laplace type operators are a special case of non-Laplace type oper- 
ators whose leading symbol is built in a universal, polynomial way, using tensor 
product and contraction from the Riemannian metric, its inverse, together with (if 
applicable) the volume form and/or the fundamental tensor-spinor. These opera- 
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tors act on sections of spin-tensor bundles. These bundles may be characterized 
as those appearing as direct summands of iterated tensor products of the tangent, 
the cotangent and the spinor bundles (see sect 2.1). Alternatively, they may be 
described abstractly as bundles associated to representations of the spin group. 
These are extremely interesting and important bundles, as they describe the fields 
in quantum field theory. The connection on the spin-tensor bundles is built in a 
canonical way from the Levi-Civita connection. The symbols of natural opera- 
tors are constructed from the jets of the Riemannian metric, the leading symbols 
being constructed just from the metric. In this case, even if the leading symbol 
is not scalar, its determinant is a polynomial in |^|^ = gi^^{x)^^^y, and, therefore, 
its eigenvalues are functions of |^| only. This allows one to use the Riemannian 
geometry and simplifies the study of such operators significantly. 

For non-Laplace operators on manifolds without boundary even the invariant 
A4 is not known, in general (for some partial results see [13, 10, 12] and the 
review [14]). For natural non-Laplace type differential operators on manifolds 
without boundary the coefficients Aq and A2 were computed in [13]. For general 
non-Laplace type operators they were computed in [10, 12]. 

The primary goal of the present work is to generalize this study to general 
non-Laplace type operators on manifolds with boundary. We introduce a "non- 
commutative" Dirac operator as a first-order elliptic partial differential operator 
such that its square is a second-order self-adjoint elliptic operator with positive 
definite leading symbol (not necessarily of Laplace type) and study the spectral 
asymptotics of these operators with Dirichlet boundary conditions. 

This paper is organized as follows. In sect. 2 we describe the construction 
of non-Laplace type operators. In sect. 2.1 we define natural non-Laplace type 
operators in the context of Stein- Weiss operators [21]. In sect. 2.2 we describe 
a class of non-Laplace type operators that appear in matrix geometry following 
[10, 12]; we develop what can be called the non-commutative exterior calculus 
and construct first-order and second-order invariant differential operators. In sect. 
2.3 we describe the general setup of the Dirichlet boundary value problem for 
such an operator and introduce necessary tools for the analysis of the ellipticity 
condition. In sect. 3 we review the spectral asymptotics of elliptic operators both 
from the heat kernel and the resolvent point of view. In sect. 4 we develop a 
formal technique for calculation of the heat kernel asymptotic expansion. In sect. 
4. 1 the interior coefficients Aq and A2 are computed (which are the same as for 
the manifolds without boundary), and in the sect. 4.2 we compute the boundary 
coefficient Ai. 
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2 Non-Laplace Type Operators 
2.1 Natural non-Laplace type Operators 

Natural non-Laplace type operators can be constructed as follows [21]. Let M 
be a smooth compact orientable n-dimensional spin manifold (with or without 
boundary). Let S be the spinor bundle over a spin manifold M and 

'V = TM(8>---(8>TM^T*M^---(8>T*M^S (2.1) 

be a spin-tensor vector bundle corresponding to a representation of the spin group 
Spin(n) and 

V : C^CV) C°°(T*M (8) 'V) (2.2) 
be a connection on "V. Then the decomposition 

rM(8)^ = nVi©---©nV, (2.3) 

of the bundle T*M (g) into its irreducible components 'Wi, 'Ws defines the 
projections 

Pj-.TM®^ ^^Vj (2.4) 
and the first-order diff'erential operators 

Gj = PjV : C"'(^) ^ C^CWj), (2.5) 

called Stein- Weiss operators (or simply the gradients). The number s of gradients 
is a representation-theoretic invariant of the bundle 'V. 

Then every first-order Spin(R) -invariant differential operator 

D : C^CV) C^CV) (2.6) 

is a direct sum of the gradients 

D = ciGi + ■■■ + CsG, = PV , (2.7) 

where cj are some real constants and 

s 

P = J]cjPj, (2.8) 



9 



and the second-order operators 

L : C~(T) C~(^) (2.9) 

defined by 

L = D*D = V*P^V = c]G)Gj (2.10) 

;=i 

are natural non-Laplace type operators. If all cj 0, then L is elliptic and has a 
positive definite leading symbol. 

2.2 Non-commutative Laplacian and Dirac Operator in Ma- 
trix Geometry 

Let M be a smooth compact orientable n-dimensional spin manifold with smooth 
boundary dM. We label the local coordinates jc^ on the manifold M by Greek 
indices which run over 1, . . . , n, and the local coordinates Jc' on the boundary dM 
by Latin indices which run over 1. We use the standard coordinate 

bases for the tangent and the cotangent bundles. The components of tensors over 
M in the coordinaate basis will be labeled by Greek indices and the components 
of tensors over dM in the coordinate basis will be labeled by Latin indices. We 
also use the standard Einstein summation convention for repeated indices. 

Let S be now an arbitrary A^-dimensional complex vector bundle over M (non 
necessarily the spinor bundle) with a positive definite Hermitean inner product 
( , ), tS* be its dual bundle and End(iS) be the bundle of linear endomorphisms 
of the vector bundle S . Further, let Aut(<S) be the group of automorphisms of the 
vector bundle S and G(S) be the group of unitary endomorphisms of the bundle 
S. We will call the unitary endomorphisms of the bundle S simply gauge trans- 
formations. 

Let TM and T*M be the tangent and the contangent bundles. We introduce 
the following notation for the vector bundles of vector- valued and endomorphism- 
valued p-forms and p-vectors 

Ap = (APT*M)®S, A" = (apTM)®S, (2.11) 

Ep = (apT*M) (8) End(<S) , = {A^TM) End(<S) . (2. 12) 

We will also consider vector bundles of densities of different weights over the 
manifold M. For each bundle we indicate the weight explicitly in the notation of 
the vector bundle; for example, S[w\ is a vector bundle of densities of weight w. 
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Since M is orientable there is the standard volume form vol = dx = dx^ A 
• • • A dx" given by the standard Lebesgue measure in a local chart. The volume 
form is, of course, a density of weight 1, and, hence, is a section of the bundle 
^^[l]. The components of the volume form in a local coordinate basis are given 
by the completely anti-symmetric Levi-Civita symbol e^^,„^^. The n-vector dual to 
the volume form is a density of weight (-1) and, hence, is a section of the bundle 
E"{-Y\. Its components are given by the contravariant Levi-Civita symbol e^^-^^". 
These objects naturally define the maps 

e : AP[w] A„_p[w + 1] , s : Ap[w] A"-p[w - 1] . (2.13) 

It is not difficult to see that 

ss = ss = (-iy^"-Phd. (2.14) 

Further, we define the diffeomorphism- invariant -inner product on the space 
C~ (Ap 5]) of smooth endomorphism- valued p-form densities of weight ^ by 



(ilf,(p)= dx {ilf{x),(p{x)) . (2.15) 



M 



The completion of in this norm defines the Hilbert space O {Ap 

Suppose we are given a map 

F : T*M End(»S) (2.16) 

determined by a self-adjoint endomorphism- valued vector 
F G C°° {TM (g) End(<S)[0]), which is described locally by a matrix-valued vector 
F''. Let us define an endomorphism- valued tensor a G C°° {TM ® TM ® End(iS)[0]) 
by 

ai^u^i) = \ immi) + r(^2)r(^i)] , (2.17) 
where ^1,^2 ^ T*M. Then a is self-adjoint and symmetric 



ai^u^i) = a(^2,^i) , a(^i,6) = • (2-18) 

One of our main assumptions about the tensor a is that it defines an isomor- 
phism 

a:TM®S^TM®S. (2.19) 
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Let us consider the endomorphism 

H{x,^) = a{^,^) = {Tm\ (2.20) 

with X e M, and ^ e T*M being a cotangent vector. Our second assumption is 
that this endomorphism is positive definite, i.e. 

H(x,^) > (2.21) 

for any point x of the manifold M and ^ ^ 0. This endomorphism is self-adjoint 
and, therefore, all its eigenvalues are real and positive for ^ ^ 0. We call the 

endomorphism- valued tensor a the non-commutative metric and the components 
of the endomorphism- valued vector Y the non-commutative Dime matrices. 
This construction determines a collection of Finsler geometries [9, 12]. As- 
sume, for simplicity, that the matrix H(x,^) = a{^,^) has distinct eigenvalues: 
h(a){x, 0,a= l,...,N. Each eigenvalue defines a Hamilton- Jacobi equation 

h^a)(x,dS) = mf^^, (2.22) 

where m^a) are some constants, a Hamiltonian system 

dxf 1 d 



(the coefficient 1/2 is introduced here for convenience) and a positive definite 
Finsler metric 



Moreover, each eigenvalue is a positive homogeneous function of ^ of degree 
2 and, therefore, the Finsler metric is a homogeneous function of ^ of degree 0. 
This leads to a number of identities, in particular, 

and 

i^=4>,^)^v. (2.27) 
Next, one defines the inverse (covariant) Finsler metrics 

g(a)/.v(^,%(:)(^,^) = 5^ (2.28) 
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the interval 

ds^a) = 8(a)Mv(x, x) dxf'dx" , (2.29) 

connections, curvatures etc (for details, see [40]). Thus, a non-Laplace type oper- 
ator generates a collection ofFinsler geometries. 
The isomorphism a naturally defines a map 

A:hp^hP, (2.30) 

by 

{A^pT"-'''' = A^'-^''^'-''''^vr-v, > (2.31) 

where 

j^t^vi^.w-y,, = glM^^ . . . gfh^gVn . . . . . . a"pPp , (2.32) 

and the square brackets denote the complete antisymmetrization over the indices 
included. We will assume that these maps are isomorphisms as well. Then the 
inverse operator 

:AP^Ap, (2.33) 

is defined by 

(A-V)....., = (A-i)^,..^^,,..,^^^--" , (2.34) 
where A~^ is determined by the equation 

(A~^) AVi-Vpai-ap cai _ _ _ cap ^2 35^ 

This can be used further to define the natural inner product on the space of 
p-forms Ap via 

<^^, cp) = 1^^,...^ A^--^.".-.^,^...,^ . (2.36) 

Let d be the exterior derivative on p-form densities of weight 

d : C°°(Ap[0]) ^ C-iAp^ilO]) (2.37) 

and d be the coderivative on p-vector densities of weight 1 

d = (-IfP^^Ms : C"'(A^[1]) ^ C^iAP-^ll]) . (2.38) 

These operators are invariant diff'erential operators defined without a Riemannian 
metric. They take the following form in local coordinates 

{d(p)^y..i,p,, = (p + l)5|^i(^^2...^^] , (2.39) 

(dcpf'-""-' = d^i/"'-""-' . (2.40) 
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and 

given locally by 



Now, let BeC°° {T*M (g) End(iS)[0]) be a smooth anti-self-adjoint endomor- 
phism- valued connection 1-fonn on the bundle S, defined by the matrix- valued 
covector !B^. Such a section naturally defines the maps: 

S : A, [i] ^ A,,i [i] (2.41) 
B = i-iyP^^sSs : AP [i] ^ A"-' [i] (2.42) 

(Scp),,..,^^, = {p+m^,ip,,..,^^,^, (2.43) 

(Sipf' "^"-' = Sf.i^^'"'^"-' . (2.44) 

Finally, we introduce a self-adjoint non-degenerate endomorphism- valued den- 
sity p e C°° (End(»S) [^]) of weight Then has weight 1 and plays the role of 
a non-commutative measure. 

This enables us to define the covariant exterior derivative of p-form densities 
of weight J 

2):C-(A4i])^C-(A,,i[i]). (2.45) 
and the covariant coderivative of p- vector densities of weight | 

£> = (-lyp+^sDs : C~ (a^ [i]) ^ C~ (a''-^ [i]) , (2.46) 

by 



D=p(d + B)p-' , (2.47) 

D =p-\d + S)p. (2.48) 

These operators transform covariantly under both the diffeomorphisms and the 
gauge transformations. 

The formal adjoint of the operator D 

l):C-(A4i])^C-(A,_i[i]), (2.49) 

has the form 

D = -A-^p-\d + S)pA , (2.50) 

By making use of these operators we define a second-order operator (that can be 
called the non-commutative Laplacian) 

A: C-[A,[W)^C-[a,[W), (2.51) 
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by 

A = -1)£)-D1). (2.52) 
In the special case p = the non-commutative Laplacian A reads 

^ = p^\d + B)pAp{d + B)p-\ (2.53) 

which in local coordinates has the form 

A = p-\d, + B,)pa>''p{d, + S,)p-i . (2.54) 

Next, notice that the endomorphism-valued vector T introduced above natu- 
rally defines the maps 

r : C~ (a'' [^]) ^ (a^^i [i]) (2.55) 

and 

f = {-ir^'sTS : (Ap [i]) ^ (a,_i [i]) (2.56) 

given locally by 

{Tipf'-^'p-' = (p + l)T^^(f^^-^'p*^^ (2.57) 

(f^)......,_. = r^^^....,^., . (2.58) 

Therefore, we can define a first-order invariant diff'erential operator (that can 
be called the non-commutative Dirac operator) 

D : C- (a, [i]) C- (A, [i]) (2.59) 

by 

D = iYD = ifpid + S)p-^ , (2.60) 
where, of course, / = V-T. The formal adjoint of this operator is 

D = iA-^DFA = iA-^p-\d + !B)pTA . (2.61) 

These operators can be used to define second order diff'erential operators DD and 
DD. 

In the case p = these operators have the following form in local coordinates 

D = /rV(5^+S>-\ (2.62) 
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D = ip-\dy + !By)pr , (2.63) 
and, therefore, the second-order operators D D and D D read 

DD = -T''p(d^+S^)p-\dy + !By)pr, (2.64) 

DD = -p-\dy + !By)prr^'p(d^ + Si,)p-K (2.65) 



In the present paper we will primarily study the second-order operators A,DD 
and D D in the case p = 0, that is, 

A,DD,DD:C^{S[^])^C'"{S[^]) . 

These are all formally self-adjoint operators by construction. This means that they 
are symmetric on smooth sections of the bundle S with compact support in the 
interior of M (that is, sections that vanish together with all their derivatives on the 
boundary dM). 

The leading symbols of all these operators are equal to the matrix H(x, ^) = 
a(^,^), i.e. 

o-l(A; X, = o-d D D; x, = o-l(DD; x, ^) = H(x, ^) = a(^, ^ , (2.66) 

where ^ e T*M. By our main assumption about the non-commuting metric the 
leading symbol is self-adjoint and positive definite in the interior of the manifold. 
Therefore, the leading symbol is invertible (or elliptic) in the interior of M. Notice 
that the leading symbol is non-scalar, in general. That is why such operators are 
called non-Laplace type operators. 

2.3 Elliptic Boundary Value Problem 

Let us consider a neighborhood of the boundary dM in M. Let x = be the 
local coordinates in this neighborhood. The boundary is a smooth hypersurface 
without boundary. Therefore, there must exist a local diffeomorphism 

r = r(x) x' = x'(x) , /=!,...,«- 1, (2.67) 

and the inverse diffeomorphism 

jt^ = jc^(r,x), ju=l,...,n, (2.68) 



such that 



r(x) = 



for any x G dM , 



(2.69) 
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r(x) > for any x ^ dM , (2.70) 

and the vector dr = d/dr is transversal (nowhere tangent) to the boundary dM. 
Then the coordinates x' are local coordinates on the boundary dM. 
Let 6 > 0. We define a disjoint decomposition of the manifold 

M = Mi„tUMb„d, (2.71) 

where 

Mbnd = {xeM\ r(x) < 6] (2.72) 
is a (^-neighborhood of the boundary and 

Mi,, = M\M^,<i (2.73) 

is the part of the interior of the manifold on a finite distance from the boundary. 

For r = 0, that is, x e dM, the vectors {di = d/dx'} are tangent to the boundary 
and give the local coordinate basis for the tangent space T^dM. The set of vectors 
{dr, di} gives the local coordinate basis for the tangent space T^M in Mbnd- Simi- 
larly, the 1 -forms dx' determine the local coordinate basis for the cotangent space 
T*dM, and the 1 -forms {dr, dx'} give the local coordinate basis for the cotangent 
space r*M in Mbnd- 

We fix the orientation of the boundary by requiring the Jacobian of this diffeo- 
morphism to be positive, in other words, for any x G Mbnd 

J(x) = vol (dr, du..., K-i) > . (2.74) 

Let 6 C°°{TM\X\) be a smooth vector density of weight 1. Then Stokes' 
Theorem has the form 

Jdxd(p = J dxN((p), (2.75) 

M dM 

where A'^ is a 1-form defined by 

N((p) = vol ((p,du..., dn-i) = jdr((p) 

dx"' dx'"-' „ I dr ^ ^^^^^ 



Notice that this formula is valid for densities, and there is no need for a Rieman- 
nian metric. 
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We will study in the present paper, for simplicity, the Dirichlet boundary con- 
ditions 

^L = 0- (2-77) 

By integration by parts it is not difficult to see that all operators A, D D and D D 
are symmetric on smooth sections of the bundle ^S^^j satisfying the boundary 
conditions. One can show that these operators are essentially self-adjoint, that is, 
their closure is self-adjoint and, hence, they have unique self-adjoint extensions to 

mil)- 

Let L be one of the operators D D,D D , A with the Dirichlet boundary con- 
ditions. Our primary interest in this paper is the study of elliptic boundary value 
problems. EUipticity means invertibility up to a compact operator in appropriate 
functional spaces (see, for example, [19, 31, 29]). This is, roughly speaking, a 
condition that implies local invertibility. For a boundary value problem it has two 
components: i) in the interior of the manifold, and ii) at the boundary. 

An operator L is elliptic in the interior of the manifold if for any interior point 
X e M and for any nonzero cotangent vector ^ e T*M, ^ 0, its leading symbol 
<jl{L', X, ^) is invertible. Since all operators DD, DD and A all have positive 
leading symbols, namely H(x, ^), they are elliptic in the interior of the manifold. 

At the boundary dM of the manifold we use the coordinates (r, x) and define a 
spfit of the cotangent bundle T*M = R © T*dM, so that ^ = (^^) = (co, |) G T;M, 
where a> e R and | = (^j) e T;dM. 

Let /I e C \ R+ be a complex number that does not lie on the positive real axis 
and H(r, x, oj, |) be the leading symbol of the operator L. We substitute r = and 
CO -idr and consider the following second-order ordinary differential equation 
on the half-line, i.e. r G R+, 

[H(0,x,-idr,i)-A]<p = 0, (2.78) 

with an asymptotic condition 

lim = . (2.79) 

r— >oo 

Let S = Sl'jM be the restriction of the vector bundle S to the boundary. The 
operator L with Dirichlet boundary conditions is elliptic if for each boundary point 
X 6 dM, each | g T^dM, each i g C \ R+, such that | and A are not both zero, 
and each f e C°° (s [\]) there is a unique solution (p{A, r, ^) to the equation (2.78) 
subject to the asymptotic condition (2.79) at infinity and the boundary condition 
atr = 

cp(A,0,i) = f. (2.80) 



18 



We have 

H(0,x,oj,i) = [A{x)to + Cix,i)f 

= A^(x)o/ +B(x,i)oj + C^{x,i), (2.81) 

where A, B, and C are self-adjoint matrices defined by 

A(x) = Y(dr) , C(x, |) = r(JxOli ' (2-82) 

B(x, I) = A(x)C(x, I) + C(x, l)A(x) . (2.83) 
Then the diff'erential equation (2.78) has the form 

(-A^dl - iBdr + - Al) ^ = . (2.84) 

We notice that the matrix |^(Aa) + C)^ - /llj is non-degenerate when co is real 
and X and | are not both zero, i.e. {A, |) (0, 0). Moreover, when /I is a nega- 
tive real number, then this matrix is self-adjoint and positive definite for real cj. 
Therefore, we can define 

oo 

/J 
^e^R,(co,l), (2.85) 

— oo 

where 

R;^(co, I) = [[A(x)oj + C(x, i)f - Aiy' . (2.86) 
The matrix ^(A,y, |) is well defined for any y e R. It: i) vanishes at infinity. 

Urn (^{A,y,l) = 0, (2.87) 

y—¥±oo 

ii) satisfies the symmetry relations 



^(A,y,^) = 0(2, -y,^) , ^(A,y, -^) = 0(A, -y,^) , (2.88) 
iii) is homogeneous, i.e. for any t > 0, 

<^^^,^^ty,^^ = t"^^iA,y,i), (2.89) 
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iv) is continuous at zero with a well defined value at y = 

oo 

-^Mco,i), (2.90) 

— oo 

v) has a discontinuous derivative dy(^(A, y,^) aty = with a finite jump. 

We also notice that the matrix Oo(/l, |) is an even function of | and is self- 
adjoint for real A, i.e. 



Oo(A, -a = Oo(^,a , Oo(^a = ^o(^,^) . (2.91) 

Moreover, for real negative A the matrix Oq is positive and, therefore, non-degenerate. 
More generally, it is non-degenerate for ReA < w, where w is a sufficiently large 
negative constant. 

In an important particular case, when B = AC + CA = 0, one can compute 
explicitly 

m,y,l) = ^A-V"'e"^'''A-' , Oo(A,|) = ^A-yU-i , (2.92) 



where /u = ^/A-\C^ - Al)A-\ defined as an analytical continuation in /I of a 
positive square root of a self-adjoint matrix when /I 6 R_. 

One can prove now that the eq. (2.84) with initial condition (2.80) and the 
asymptotic condition at infinity (2.79) has a unique solution given by 

ip(A, r, I) = 0(^, r, |)[Oo(^, hv'f . (2.93) 

Thus, the Dirichlet boundary value problem for our operator is elliptic. 



3 Spectral Asymptotics 
3.1 Heat Kernel 

Let L be a self-adjoint elliptic second-order partial diff'erential operator acting on 
smooth sections of the bundle over a compact manifold M with boundary 
dM with positive definite leading symbol and with some boundary conditions 



Bip\aM = , 



(3.1) 
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with some boundary operator B. It is well known that such an operator has a 
discrete real spectrum {/Ifcl^j bounded from below [29], i.e., 

< I2 < • • • < 4 < 4+1 < • • • . (3.2) 

Furthermore: i) each eigenspace is finite-dimensional, ii) the eigensections are 
smooth sections of the bundle ^j, and iii) the set of eigensections {(Pk}'k=i forms 
an orthonormal basis in (^S 
For ? > the heat semigroup 

exp(-?L) : [S [i]) ^ [S [i]) (3.3) 

is a bounded operator. The integral kernel of this operator, called the heat kernel, 
is given by 

00 

U{t, X, x') = 2 e-'^'ifik <» , (3.4) 
k=i 

where each eigenvalue is counted with its multiplicity. The heat kernel satisfies 
the heat equation 

(dt + L)U(t;x,x') = (3.5) 

with the initial condition 

U{0^;x,x') = 6{x,x'), (3.6) 
where 5(x, x') is the Dirac distribution, as well as the boundary conditions 

BMt;x,x') =0, (3.7) 

x€dM 

and the self-adjointness condition 



U{t;x,x') = U{t;x',x). (3.8) 

The heat kernel U{t) = exp(-rL) is intimately related to the resolvent G(A) = 
(L - Ay^. Let /I be a complex number with a sufficiently large negative real part, 
Re A « 0. Then the resolvent and the heat kernel are related by the Laplace 
transform 

CO 

G(A) = J dte'^Uit), (3.9) 
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W+l(x. 

U(t) = J- f dAe-'-'GiA), (3.10) 
2m J 

w—icc 

where w is a sufficiently large negative real number, w « 0. 
The resolvent satisfies the equation 

(L - AI)G(A; X, x') = 6{x, x') (3.11) 

with the boundary condition 

B,G{A-x,x%^^^ = Q, (3.12) 
and the self-adjointness condition 



G{A;x,x') = G{A;x',x). (3.13) 
The integral kernel of the resolvent reads 

G{A- X, x') = yrri ® ' (^-^"^^ 



k=i ^^-^ 



where each eigenvalue is counted with its multiplicity. 



For ? > the 
product bundle S 



leat kernel lJ{t\ x, x') is a smooth section of the exterior tensor 
El S* [^j; that is, it is a two-point density of weight ^ at each 
point. In particular, it is a smooth function near the diagonal of M x M and has a 
well defined diagonal value U (t; x, x). The diagonal is, of course, a smooth section 
of the bundle [1], a density of weight 1. 

Moreover, the heat semigroup is a trace-class operator with a well defined 
-trace 

Tr i2 exp(-tL) = ^ dx tr sU(t; x,x) , (3.15) 

M 

where tr 5 is the trace over the fiber vector space S of the vector bundle S. The 
trace of the heat kernel is a spectral invariant of the operator L since 

00 

Tr L2 exp(-tL) = ^ e''^' . (3.16) 



Since the diagonal is a density of weight 1 the trace Tr 12 exp(-tL) is invariant 
under difi'eomorphisms. 
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This enables one to define other spectral functions by integral transforms of 
the trace of the heat kernel. In particular, the zeta function, ^(L; s. A), is defined as 
follows. Let /I be a complex parameter with Re /I < /li , so that the operator (L - A) 
is positive. Then for any s e C such that Re* > n/2 the trace of the operator 
(L - A)~^ is well defined and determines the zeta function. 



oo 



^(L; s. A) = Tr l2(L - A)'' = — J dt f e'^ Tr ji exp(-?L) . (3.17) 



The zeta function enables one to define further the regularized determinant of the 
operator (L - X) by 

—^(L;s,A) = -logDet(L-^). (3.18) 

OS s=0 

There is an asymptotic expansion as ? ^ of the trace of the heat kernel [29] 
(for a review, see also [2, 5, 6, 7, 43]) 

oo 

Tvl2 exp(-?L) ~ (Anr^^ ^ ^'^''"^^^MQ . (3.19) 

k=Q 

The coefficients Ak(L), called the global heat invariants, are spectral invariants 
determined by the integrals over the manifold M and the boundary dM of some 
scalar densities ak(L; x) and bk(L; x), called local heat invariants, viz. 

Ak(L) = J dx ak(L; x) + J dx bk(L; x) . (3.20) 

M dM 

The local heat invariants ak{L; x) and bk{L; x) are constructed polynomially from 
the jets of the symbol of the operator L; the boundary coefficients bk depend, of 
course, on the boundary conditions and the geometry of the boundary as well. 

Contrary to the heat kernel, the resolvent is singular at the diagonal and does 
not have a well defined trace. However, the derivatives of the resolvent do. Let 
m > njl. Then the trace i2{d^)'"G{A) is well defined and has the asymptotic 
expansion as /I ^ -oo 

ffn °° 

Ti-i^^^^(^) ~ (4^)""^' V r [(/: - n + 2m + 2)/2] {-At-^-^"'-^^I^Ak{L) . 

(3.21) 

Therefore, one can use either the heat kernel or the resolvent to compute the coef- 
ficients Ak. 
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3.2 Index of Noncommutative Dirac Operator 

Notice that the operator A can have a finite number of negative eigenvalues, whereas 
the spectrum of the operators D D and D D is non-negative. Moreover, one can 
easily show that all non-zero eigenvalues of the operators D D and DD are equal 

Ak( DD) = Ak(D D) if Xk{DD)>Q. (3.22) 

Therefore, there is a well defined index 

Ind(D) = dim Ker {DD)- dim Ker {D D) , (3.23) 

which is equal to the diff'erence of the number of zero modes of the operators D D 
and DD. 

This leads to the fact that the difi'erence of the heat traces for the operators D D 
and D D determines the index 

Tr L2 exp(-t DD)-Ty ^2 exp(-? DD) = Ind(D) . (3.24) 

This means that the spectral invariants of the operators D D and DD art equal 
except for the invariant A„ which determines the index 

Ak(DD)=Ak(DD) for ki^n, (3.25) 

and 

A„( DD)- A„(D D ) = (4nr'Hnd(D) . (3.26) 

Thus, for n > 2 the spectral invariants Aq, Ai and of the operators D D and 
D D are equal. Therefore, we can pick any of these operators to compute these 
invariants. Of course, the spectral invariants of the noncommutaative Laplacian A 
are, in general, diff'erent. However, since the operators D D and D D have the same 
leading symbol as the operator A there must exist a corresponding Lichnerowicz- 
Weitzenbock formula (for the spinor bundle see, for example, [18]), which means 
that the spectral invariants of these operators must be related. 

4 Heat Invariants 
4.1 Interior Coefficients 

The heat kernel in the interior part is constructed as follows. We fix a point xq e 
Mint in the interior of the manifold and consider a neighborhood of xq disjoint 
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from the boundary layer Mbnd covered by a single patch of local coordinates. We 
introduce a scaling parameter e > and scale the variables according to 

jt^ ^ jc^ + - jc^) , x''^ ^ jc^ + - jc^) , t ^ eh , (4.1) 

so that 

d^^^-d^, d,^\d,. (4.2) 
Then the differential operator L(x, d) is scaled according to 

oo 

L^L,~YjE'-'i};:\ (4.3) 

k=0 

where L™^ are second-order differential operators with homogeneous symbols. 
Next, we expand the scaled heat kernel in Mint, which we denote by i7g"' in a 
power series in s 

oo 

k=0 

and substitute into the scaled version of the heat equation. By equating the like 
powers of e we get an infinite set of recursive differential equations determining 
all the coefficients Up. 

The leading order operator is an operator with constant coefficients deter- 
mined by the leading symbol 

l}^' = H(xo,-id). (4.5) 
The leading-order heat kernel U^^ can be easily obtained by the Fourier transform 

Up(t;x,x^)= f ^g^(-^W^(-o,f) ^ (4.6) 
J (27r)" 

R" 

where ^(x - x') = ^^(x^' - x'^'). 

The higher-order coefficients Up, k > 1, are determined from the recursive 
equations 

k 

(d, + Lf)Up = -Y, l^Utj > (4.7) 

7=1 

with the initial condition 

Up{Q)-x,x') = Q>. (4.8) 
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This expansion is nothing but the decomposition of the heat kernel into the 
homogeneous parts with respect to the variables (x - xq), (x' - xq), and V^- That 
is, 

[/f (r; X, X') = UT (l ; xo + ;co + ^^^^ j . (4.9) 

In particular, the heat kernel diagonal at the point xq scales by 

Ut\t; xo, xo) = fi-"^'^ [/^' (1 ; xo, xo) . (4. 10) 

To compute the contribution of these coefficients to the trace of the heat kernel 
we need to compute the integral of the diagonal of the heat kernel U"^^{t; x, x) over 
the interior part of the manifold Mint. By using the homogeneity property (4.10) 
we obtain 

r dxiVsU'''\t;x,x) ~ ^(^-«)/2 f Jjctr^f/f (l;x,x) . (4.11) 

^ k-O 

Next, we take the limit as 6 ^ 0. Then the integrals over the interior part Mint be- 
come the integrals over the whole manifold M and give all the interior coefficients 
ak(L) in the global heat kernel coefficients Ak(L). 

Instead of this rigorous procedure, we present below a pragmatic formal ap- 
proach that enables one to compute all interior coefficients in a much easier and 
compact form. Of course, both approaches are equivalent and give the same an- 
swers. 

First, we present the heat kernel diagonal for the operator L = D D in the form 
ir^\t;x,x)= [ -^e-^-exp(-?DD)e^\ (4.12) 

where ^x = ^/^^f, which can be transformed to 

U^\t,x,x) = J ^cxp[-t(H + K + DD)]-1, (4.13) 

w 

where H = [r(^)]^ is the leading symbol of the operator D D, and is a first-order 
self -adjoint operator defined by 



K = -mD- Dm. 



(4.14) 



26 



Here the operators in the exponent act on the unity matrix I from the left. 
By changing the integration variable ^ —> t~^^^^ we obtain 

U"^\t; X, x) = (Ant)-"'^ J exp (-H - yft K - tD d) -l. (4. 15) 

R" 

Now, the coefficients of the asymptotic expansion of this integral in powers of t^l^ 
as ? ^ determine the interior heat kernel coefficients ak{L) via 



tr s [/'"'(f; X, x) ~ (An)-"'^ ^ t^^~"^'^ak{L) 
By using the Volterra series 

exp(A + 5) = + 2 j Jr^t I dTk-i--- I dT 



(4.16) 



k=0 



1 n 



X 



we get 







1 



1 T2 



1 



(4.17) 



-/ 



dTie 



-(l-ri)ff 



DDe 



-TlH 



+ 0(t'). 



(4.18) 



Now, since K is linear in ^ the term proportional to t^^^ vanishes after integra- 
tion over ^. Thus, we obtain the first three interior coefficients of the asymptotic 
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expansion of the heat kernel diagonal in the form 

«o(^) = ^%''se-\ (4.19) 

R" 

ai(L) = 0, (4.20) 

1 rz 

R" 

1 







(4.21) 



4.2 Boundary Coefficients 

On manifolds with boundary, as far as we know, the coefficients Ak have not been 
studied at all, so, even Ai is not known. In the present paper we are going to com- 
pute the coefficient A 1 on manifolds with boundary for the operators DD and DD. 
The coefficient Aq is, of course, the same as for the manifolds without boundary. 
We will follow the general framework for computation of the heat kernel asymp- 
totics outlined in [1 1, 16]. 

The procedures for the resolvent and the heat kernel are very similar. One can, 
of course, use either of them. We will describe below the construction of the heat 
kernel. 

The main idea can be described as follows. Recall that we decomposed the 
manifold into a neighborhood of the boundary Mbnd and the interior part Mi^. 
We can use now different approximations for the heat kernel in different domains. 
Strictly speaking one has to use 'smooth characteristic functions' of those domains 
(partition of unity) to glue them together in a smooth way. Then, one has to control 
the order of the remainder terms in the limit t ^ and their dependence on 6 
(the size of the boundary layer). However, since we are only interested in the trace 
of the heat kernel, this is not needed here and we will not worry about such subtle 
details. We can compute the asymptotic expansion as ? ^ of the corresponding 
integrals in each domain and then take the limit (5^0. 

The origin of the boundary terms in the heat trace asymptotics can be ex- 
plained as follows. The heat kernel of an elliptic boundary value problem in Mbnd 
has exponentially small terms like exp(-r^/0 as ? ^ 0. These terms do not con- 
tribute to the asymptotic expansion of the diagonal of the heat kernel as ? ^ 0. 
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However, they behave like distributions near the boundary (recall that r > in- 
side the manifold and r = on the boundary). Therefore, the integral over Mbnd, 

6 

more precisely, the limit lim^^o / dx J dr(. . . ) does contribute to the asymptotic 

dM 

expansion of the trace of the heat kernel with coefficients in form of integrals over 
the boundary. It is this phenomenon that leads to the boundary terms in the global 

heat invariants. 

The heat kernel in the boundary layer Mbnd is constructed as follows. We fix 
a point xq e dM on the boundary and choose coordinates as described above 
in section 2.2. Let £ > be a positive real parameter. We use it as a scaling 
parameter; at the very end of the calculation it will be set to 1. Now we scale the 
coordinates according to 

x^ ^ xi + s(x^ - x^) , x'j ^xl + s(x'^ - x^) , (4.22) 

r ^ sr , r ^ sr , t ^ s^t . (4.23) 
The difi'erential operators are scaled correspondingly by 

dj ^ -dj , dr ^ -dr , dt ^ ^d, . (4.24) 

Let L(r,x,dr,d) be the operator under consideration. The scaled operator, 
which we denoted by Lg, has the following formal power series expansion in s 

oo 

Lh^L,~2e^-24"^ (4.25) 

*:=0 

where are second-order differential operators with homogeneous symbols. The 
leading order operator is determined by the leading symbol 

l!^"'' = H(0, xo, -idr, -id) . (4.26) 

This is a diff"erential operator with constant coefficients. 

Next, we expand the scaled heat kernel in Mbnd, which we denote by i/^"'^ in a 
power series in s 

oo 

k=Q 

and substitute into the scaled version of the heat equation and the boundary con- 
ditions. By equating the like powers of e we get an infinite set of recursive difi'er- 
ential equations determining all the coefficients C/]^"'^. 
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The leading-order heat kernel U^'^^ is determined by the equation 

(5f + L^"^)f/^"'' = (4.28) 

with the initial condition 

U^^^O; r, X, /, x) = 6{r - r')6{x, x) , (4.29) 
the boundary condition 

Ul"^{t; 0, X, r', x') = U^''\t; r, x, 0, x') = , (4.30) 
and the asymptotic condition 

lim U^"\t; r, x, /, x) = lim U^''\t; r, x, r',x) = 0. (4.3 1) 



The higher-order coefficients C/^"^, ^ > 1, are determined from the recursive 
equations 

k 

(d, + L^^^Wl^" = -J] LfU^fj , (4.32) 
with the initial condition 

Ul''\0;r,x,r',x') = 0, (4.33) 

the boundary condition 

U^^"\t; 0, X, r', x') = U''^\t\ r, x, 0, x) = 0, (4.34) 
and the asymptotic condition 

lim [/^"'^(r; r, x, /, x) = lim U^''\t; r, x, /, x) = 0. (4.35) 



This expansion is nothing but the decomposition of the heat kernel into the 
homogeneous parts with respect to the variables (Jc - xq), (x' - xq), r, f and V^- 
That is, 



(4.36) 
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In particular, the heat kernel diagonal at the point (r, xq) scales by 




(4.37) 



To compute the contribution of these coefficients to the trace of the heat kernel 
we need to compute the integral of the diagonal of the heat kernel C/'^°'*(f; r, x, r, x) 
over the boundary layer Mbnd- This heat kernel diagonal can be decomposed as 
the sum of two terms, the first coming from the standard interior heat kernel on 
manifolds without boundary (that does not satisfy the boundary conditions) and 
the second 'compensating' part, which is the crucial boundary part and whose role 
is to make the heat kernel to satisfy the boundary conditions (for more details see 
[11]). The integral of the 'boundary' part over the boundary layer in the limit when 
the size of the boundary layer goes to zero produces the boundary contributions 
to the global heat kernel coefficients A^^(L). 

By using the homogeneity property (4.37) we obtain 



where u = r/ y/t. Notice the appearance of the extra power of V? in the asymptotic 
expansion. Of course, if one takes the limit lim^^o for a finite t, then all these 
integrals vanish. However, if one takes the limit lim,_>o first for a finite 6, and then 
the limit lim^^o. then one gets finite answers for the boundary coefficients b^iL). 

4.2.1 Leading-Order Heat Kernel 

To compute the coefficient Ai we just need the leading-order heat kernel U^'^^. We 
will, in fact, be working in the tangent space R+xTj^^dM at a point Xq on the bound- 
ary and reduce our problem to a problem on the half -line. The operator Lq""* acts on 



6 






(4.38) 
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square integrable sections of the vector bundle S[j] in a neighborhood of the point 
xq. We extend the operator appropriately to the space L^(<S[^], R+, , dr dx) so 
that it coincides with the initial operator in the neighborhood of the point xq. When 
computing the trace below we set xq = x = x'. 

By using the Laplace transform in the variable t and the Fourier transform in 
the boundary coordinates x 

w+ioo ^ 

U^'^'it; r, X, r', x') = ^. ^ dX j ^^^-'^-^C-^) F(A, , r, r', |) , (4.39) 

W-ioo Rn-l 

we obtain an ordinary differential equation 

{-A^dl - iBdr + C^-M) F(A, r, r' , |) = I6{r - r') (4.40) 

where the matrices A, B and C are defined in (2.82), (2.83), and are frozen at the 
point Xq (they are constant for the purpose of this calculation), with the boundary 
condition 

F{X,Q,r',l) = F{X,r,Q,l) = Q (4.41) 

the asymptotic condition 

lim F{A, r, r', |) = lim F{A, r, r', |) = , (4.42) 

r— >oo r'— >oo 

and the self-adjointness condition 

F{A,r,r',l) = F(l, r', r, |) . (4.43) 
It is easy to see that F is a homogeneous function 




(4.44) 



We decompose the Green function in two parts, 

F = F^ + Fb, (4.45) 

where Foo is the part that is valid for the whole real line and Fg is the compensating 
term. The part Foo can be easily obtained by the Fourier transform; it has the form 

FM r, r, I) = 0(A, r - r' , |) , (4.46) 
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where 0(/l, r,|) is defined in (2.85). It is not smooth at the diagonal r = r' and 
is responsible for the appearance of the delta-function 6(r - f) on the right-hand 
side of the eq. (4.40). 

The corresponding part of the leading heat kernel is then easily computed to 

be 

UXt{t-x,^) = J _^e^(-^W^(-o,f) ^ (4.47) 

where xq = (0, xq). This part does not contribute to the asymptotics of the trace of 
the heat kernel in the limit 6 —> O.By rescaling ^ ^/ V? we obtain 

J dxtrs Ul'i(t;x,x) = {Antr'^ J dx J ^tv se-"^''^^ , (4.48) 

Mbnd Mbnd R" 

and in the limit 6^0 this integral vanishes. 

However, Fc<, does not satisfy the boundary conditions. The role of the bound- 
ary part, Fb, is exactly to guarantee that F satisfies the boundary conditions. The 
function is smooth at the diagonal r = It can be presented in the following 
form 

Fb(A, r, r', I) = -0(A, r, |)[Oo(A, hV^i^, -r' , |) . (4.49) 
4.2.2 The Coefficient Ai 

The coefficient Ai is a pure boundary coefficient that is computed by integrating 
the boundary part of the heat kernel. We have 

^ dxiisUl'^l{t;x,x) (4.50) 



^ W+lOO 



dM Rn-i 

Now, by rescaling the variables 



X^-, r^yftr, | ^ (4.51) 
t it 
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and using the homogeneity property (4.44) we obtain 

J dxtrsU^'^^(t,x,x) (4.52) 



*^bnd 

6/ ift w+ioo 



dM R"-! w-i'oo 

Therefore, the coefficient Ai is given by 

oo w+ioo 



A,=2^Jdx JJ^,Jdr J ^_e-hrsFs(A,r,r,l). (4.53) 

dM Rn-l W-/00 

Thus, finally, by using eq. (4.49), eliminating the odd functions of | (since the in- 
tegrals of them vanish), using the property (2.88) of the function O and extending 
the integration over r from -oo to +oo (since the integrand is an even function) we 
obtain 

^1 = Jdx J Ji^^iil) (4.54) 



-l)/2 

dM R"-l 



where 



oo W+/00 



^i(l) = -^Jdr J ^_e-'trs[Oo(A,i)r' (4.55) 



— oo w—loo 



x[0(A, r, |)0(^, -r, I) + 0(^, -r, |)0(^, r, |)) . 

Recall that w is a sufficiently large negative constant. 

Now, using eq. (2.90) and integrating over r we obtain finally 

w+ioo 

— e-^tr5[Oo(A,|)]-i— Oo(A,|) 

IV- (OO 

W+/00 

— e--* — logdet[Oo(^,|)]. (4.56) 

W-!00 

Thus, the problem is now reduced to the computation of the integral over A. 
This is not at all trivial because of the presence of two non-commuting matrices. 
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essentially, A'^AC + CA)A-^ and A'^C^ - AI)A-\ where the matrices A = r(x) 
and C = Y-'(x)^j are defined by (2.82). We will report on this problem in a future 
work. Here let us just mention that in the particular case when B = AC + CA = 
(for example, this is so in the case of the original Dirac operator) we get 

O -I 

tr5 [Oo(^,|)]-'^Oo(^,|) = -trs (C^ - AI)-' , (4.57) 

OA 2 

and, therefore, one can compute the integral over A to obtain 

Of course, for Laplace type operators, when [C(i,|)]^ = I^'^(i)|;|^, the integral 
can be computed explicitly, which gives the induced Riemannian volume of the 
boundary, Ai = -( s/n/l) N Yol(dM), and coincides with the standard result for 
Dirichlet Laplacian [29]. 
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